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RESUMEN 

La distribution de energia espectral de los cuasares presenta un empinamicnto abrupto del continuo a longitudes 
mas cortas que ~ 1100 A. Presentamos un nucvo modelo de extincion por polvo constituido de granos dc 
carbono cristalino. Analizamos las propiedades en extincion poco usuales de este polvo y demostramos como el 
empinamicnto abrupto obscrvado en cuasares puede ser modelizado cxitosamente por polvo hecho de carbono 
cristalino. 

ABSTRACT 

The spectral energy distribution of quasars shows a sharp steepening of the continuum shortward of ~ 1100 A. 
We present a new dust extinction model consisting of crystalline carbon grains. We analyze the unusual 
extinction properties of this dust and proceed to show how the observed UV spectral break can be successfully 
modelled by crystalline carbon dust. 
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1. INTRODUCTION 

The ultraviolet energy distribution of quasars 
is characterized by the so-called "big blue bump", 
which peaks in vF v at approximately 1000 A. 
The quasar 'composite' spectral energy distribution 
(SED) of Telfer etal. (2002, hereafter TZ02) ob- 
tained by co-adding 332 HST-FOS archived spectra 
of 184 quasars between redshifts 0.33 and 3.6, ex- 
hibits a steepening of the continuum at ~ 1100 A. A 
fit of this composite SED using a broken powerlaw 
reveals that the powerlaw index changes from —0.69 
in the ncar-UV to -1.76 in the far-UV (see TZ02). 
We label this observed sharp steepening the 'far-UV 
break'. In these proceedings, we first give a short 
overview of our classification of the HST-FOS quasar 
spectra based on the appearance of the far-UV break. 
We then describe our dust absorption model based 
on crystalline carbon grains. We illustrate the ex- 
tinction properties of this dust and proceed to com- 
pare with two quasar spectra. In a previous proceed- 
ing (Binette et al. 2005b), we described how we came 
to consider the possibility of carbon crystallite dust. 
The argumentation in support of the dust absorp- 
tion interpretation of the UV break has been fully 
described in Binette etal. (2005a, hereafter BM05). 
Complementary information about crystalline dust 
can be found in this and other proceedings (Binette 
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etal. 2005b, 2006). 

2. CLASSIFICATION OF THE FAR-UV BREAK 
APPEARANCE 

The concept of a 'composite' Spectral Energy 
Distribution (SED) that extends into the far-UV 
rests on the rcdshift effect, which allows us to ob- 
serve the far-UV region of high-zquasars, beyond 
the Lyman limit, even though our Galaxy is highly 
opaque in this wavelength region (observer- frame) . 
For the concept of a 'composite' SED obtained by 
averaging to work, the environment of quasars must 
be highly transparent beyond the Lyman limit, a 
condition which is generally satisfied for quasars at 
z q < 3.6. In what follows, we will consider the sample 
of HST-FOS spectra that H. Telfer kindly lent to us. 
These have already been corrected for Galactic dust 
extinction, for the presence of the Lyman valley due 
to unresolved Lya forest lines and for the presence 
of Lyman limit systems. Of all the spectra avail- 
able, we consider only those 106 quasars for which 
more than one grating was used since otherwise the 
Xrest coverage turned out to be too narrow to provide 
meaningful constraints to the models. Furthermore, 
only 61 of these spectra extend sufficiently in the far- 
UV (900 A) to reveal enough of the far-UV break to 
be of use to us. 

After reviewing the HST-FOS sample, some pat- 
terns emerge among the selected 61 multigrating 
spectra. The majority of them (44 quasars) resem- 
ble the composite SED of TZ02. We classify those 
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Fig. 1. Panel a (left): cartoon that illustrates in Fx the three main classes of far-UV break (continuous lines) observed 
among the 61 multigrating HST-FOS spectra that extend down to at least 900 A. The dotted line shows typical variations 
within a given class. Panel b (right): same descriptive cartoon using instead the vF v representation. 



as class (A). Others (6 quasars) show a sharper and 
much more pronounced break. We label these class 
(B). Eight other quasars present a somewhat shal- 
lower break, which appears to start further in the 
near-UV or to present a SED which is flat through- 
out the Fx region covered. A cartoon of theses three 
classes is presented in Fig. la. The dotted lines show 
typical variations within each class. Finally, only 3 
quasars did not fit any meaningful pattern (class D). 
In Fig. lb, we show the same cartoon, but in vF v 
(= AF\). We found the traditional Fx representa- 
tion of Fib. la. to be more appropriate and useful 
when fitting the UV break region. 

3. THE EXTINCTION CURVE PRODUCED BY 
CRYSTALLINE CARBON 

At the onset of our project, it became clear that 
ISM or even SMC-like dust could not reproduce a 
break as sharp as the one observed in the far-UV 
of quasars. Such a conclusion had been reached 
by Shang etal. (2005). A bibliographical research 
indicated that crystalline carbon might possess the 
required optical properties (Mutschkc etal. 2004). 
Assuming spherical grains and the Mic theory, we 
used the subroutine BHMIE (Bohrcn and Huffman 
1983) to compute the extinction curves, adopting 
the complex refraction indices n + ik as tabulated 
in Mutschke et al. (2004) and in Edwards & Philipp 
(1985) . Crystalline carbon can appear either in a 
pure form or with surface impurities. In the former 
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Fig. 2. Extinction cross section for nanodiamonds with 
and without surface impurities (Al and Dl, respec- 
tively). The small grain size regime was assumed. The 
complex refraction indices n + ik used in the derivation 
of the extinction curves are from Mutschke et al. (2004) 
for the Allende meteorite nanodiamonds and from Ed- 
wards & Philipp (1985) for the cubic (terrestrial) nan- 
odiamonds. 

case, they are referred to as terrestrial cubic dia- 
monds. In the latter case, nanodiamonds of the type 



DUST AND THE SPECTRAL ENERGY DISTRIBUTION OF QUASARS 



3 



that are found in abundance in primitive meteorites 
will represent the case of crystalline C with impu- 
rities. Convincing evidence of nanodiamonds being 
detected in emission is provided by the 3.43 and 3.53 
/im emission bands observed in two Ac/Be stars and 
one post-AGB star (see Van Kerckhoven, Tielens & 
Waelkens 2002). These bands result from the ex- 
citation of surface CH stretching modes (implying 
the existence of nanodiamonds with H surface im- 
purities). We present in Fig. [5] the extinction curves 
Al and Dl resulting from nanodiamonds with and 
without surface impurities. We assumed a powerlaw 
grain size distribution within the range 2 to 25 A. It 
turns out that within this small grain size regime, us- 
ing a log-normal distribution or altering slightly the 
above size limits do not produce any differences in 
the derived extinction curves, as explained in BM05. 

It is apparent from Fig. [2] that the steepness of 
the near- f/Uextinction rise ought to produce a sharp 
absorption break, a desired feature for the any dust 
grain model that aims at fitting the UV quasar 
break. The extinction curves in Fig. [2] were normal- 
ized in such a way that they represent the case of 
having all the carbon in the dust (assuming a solar 
C/H abundance ratio). This normalization would 
need to be scaled according to the actual but un- 
known dust-to-gas ratio appropriate to the quasar 
environment. 

4. POWERLAWS ABSORBED BY 
NANODIAMOND DUST 
Where might the putative crystalline dust be lo- 
cated? In BM05, we present calculations in which 
the dust is either intrinsic to the quasars or inter- 
galactic. In the end, we concluded that only the in- 
trinsic dust hypothesis was satisfactory. The calcu- 
lations presented here will assume the intrinsic dust 
case. 

If one fits with a powerlaw the near-UV spec- 
tra of class (A) spectra, one finds that the average 
spectral index is cxnuv = —0.44 with a dispersion 
of 0.21. The corresponding histogram is shown in 
Fig. 4. This average considers only the 21 objects, 
for which a reliable estimate of oinuv could be de- 
termined directly from the HST spectra. It is signifi- 
cantly harder than the mean value of —0.69 reported 
by TZ02 and the median value of -0.83 ± 0.04 for 
local AGN reported by Scott etal. (2004), presum- 
ably, because the softer class (C) spectra are not in- 
cluded in our average. In Fig. 5a and 5b, we present 
a powerlaw of index a v = —0.4 absorbed by nanodi- 
amond dust of Al and Dl types, respectively. The 
dust screen takes on four different thicknesses cor- 
responding to H columns (in units of 10 20 cm~ 2 ) of 
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Fig. 3. Absorption models assuming a powerlaw SED of 
index a u = —0.4 (dotted line). The nanodiamond dust 
screen corresponds to an hydrogen column of 10 20 cm~ 2 . 
The absorbed SED represented by the continuous line 
corresponds to pure cubic diamonds (Dl) while the long 
dashed line corresponds to nanodiamonds from the Al- 
lende meteorite (Al). The grey lines correspond to a 
mixture of the two flavors: dashed line: 60% Dl + 40% 
Al, dot-dashed line: 30% Dl + 70% Al. Labelled point- 
ers indicate where emission lines are expected in the rest- 
frame quasar spectra. 




Fig. 4. Distribution of the near-UV powerlaw indices 
found among the 21 class (A) quasars for which a reliable 
estimate of olnuv could be determined directly from the 
HST spectra. 

V 20 = 0.5, 1.0, 2.0 and 4.0. When V 20 « 1.0, the 




Fig. 5. Absorbed powerlaw SED assuming a constant index ot v = —0.4 (dotted line) and the following H columns in 
units of 10 20 cm -2 : N%o = 0.5, 1.0, 2.0 and 4.0. Panel a (left): terrestrial nanodiamond dust (Dl). Panel b (right): 
meteoritic nanodiamonds (Al). (fm is the fraction of Dl dust). 



absorption results in a near horizontal Fx spectral 
segment shortward of the break. The position of the 
break occurs at different wavelengths, depending on 
the dust composition (compare Fig. 5a and 5b). It 
occurs longward and shortward of Lya in the case 
of the Al and Dl dust, respectively. Furthermore, 
terrestrial diamond dust can result in a conspicuous 
absorption dip near 650 A (continuous line in Fig. 5b 
and 3), which may have some relevance to the nar- 
row absorption dip reported by Scott et al. in their 
composite SED, which was interpreted as blucshiftcd 
Nevni absorption. 

5. COMPARISON OF OUR MODELS WITH 
TWO QUASAR SPECTRA 

In Fig. 6a and 6b, we present our absorbed pow- 
erlaw fit to the class (A) quasar PG 1008+1319 and 
to the class (B) quasar Pks 1229—02, respectively. 
In the case of PG 1008+1319, it is interesting to 
notice that the UV break is followed by a far-UV 
rise, which is well modeled by our nanodiamond dust 
model, using ^20=1-2 and a dust composition con- 
sisting of 70% of Allcndc nanodiamonds and 30% 
of cubic nanodiamonds (/rji=0.30). The assumed 
spectral index is a^uv = +0.13, a value determined 
by Neugebauer etal. (1987). In the case of the class 
(B) quasar Pks 1229—02, the model shown in Fig. 6b 
assumes a^uv = —0.4 and is characterized by a col- 
umn ^20=3.6 and dust composed mostly of cubic 
diamonds (Dl) with / D i=0.90. 



It appears that the six class (B) quasars in our 
sample are not only more absorbed than class (A) 
but are extincted by dust dominated by cubic di- 
amond grains (Dl). On the other hand, the UV 
break in class (A) spectra generally requires a mix- 
ture of the two dust flavors. The absorption column 
is lesser in class (A) spectra. The break in 39 class 
(A) quasars is well reproduced using a column in the 
range 0.6 < N 20 < 1-4. 
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CyT grant 40096-F. We thank Randal Tclfer for 
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Starkmeth helped us with proofreading. 
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Fig. 6. Left panel a: rest-frame spectrum of class (A) quasar PG 1008+1319. Notice the far-UV rise shortward of the 
UV break. The model (grey line) assumes a„ = +0.13, iV2o = L2 and /m=0.3. Right panel b: rest-frame spectrum of 
class (B) quasar Pks 1229— 02. The model (grey line) assumes a v = —0.4, Y2o=3.6 and /m=0.90. This is the largest 
column encountered among classes (A+B). 
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